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Abstract 



I study bottom quark fragmentation in the Standard Model Higgs decay 
H bb, within the framework of perturbative fragmentation functions. 
I resum large collinear logarithms ln(m|^/m^) in the next-to-leading 
logarithmic (NLL) approximation, using the DGLAP evolution equations. 
Soft contributions to the MS coefficient function and to the initial condition 
of the perturbative fragmentation function are resummed to NLL accuracy 
as well. The implementation of collinear and soft resummation has a relevant 
impact on the b energy spectrum, which exhibits a milder dependence 
on factorization and renormalization scales and on the Higgs mass. I 
present some predictions on the energy distribution of 6-flavoured hadrons 
in Higgs decay, making use of data from LEP and SLD experiments to fit 
a few hadronization models. I also compare the phenomenological results 
yielded by a few processes recently provided with NLL collinear and soft 
resummations. 



1 Introduction 



The Higgs boson plays a crucial role in the Standard Model (SM) of electroweak inter- 
actions as it is responsible for the mechanism of mass generation. However, this particle 
has not yet been experimentally discovered. 

Searches for the Standard Model Higgs boson have been performed at the LEP 
collider, are currently under way at the Tevatron, and will be ultimately one of the 
main goals of experiments at the LHC. (see, for a review, e.g. [1]). In detail, the LEP 
experiments have set a lower bound on the Higgs mass at itlh > 114.4 GeV [2], mainly 
using the production channel e+e~ HZ. The Tevatron will be able to exclude a 
Higgs boson with mass lower than 130 GeV within three standard deviations [3] . Future 
experiments at the LHC will be capable of going beyond and exploring the Higgs mass 
spectrum from 100 GeV to about 1 TeV [4] . 

In order to accurately perform such searches, the use of precise QCD calculations 
will be fundamental. In this paper, I consider the decay of the Standard Model Higgs 
boson into bb pairs, i.e. H bb. In fact, the favourite discovery channel of the Higgs 
at the Tevatron consists of processes where H is produced in association with a vector 
boson, i.e. pp — > VH, where V is a Z or a W, followed by the decays H bb and 
V £ii2, ii and £2 being leptons. At the LHC, the process gg{qq) H ^ bb will 
be affected by large QCD backgrounds, which make the detection of this decay channel 
more cumbersome. However, the process H ^ bb will still play a role, in particular for 
nT'H ^ 135 GeV and Higgs production in association with tt pairs, i.e. pp — > ttH [5], with 
a W boson [6], in vector boson fusion [7]. 

Hereafter, I shall address the issue of multiple gluon radiation in H ^ bb pro- 
cesses. While fixed-order calculations are reliable enough to predict total cross sections 
or widths, differential distributions present terms, corresponding to collincar- or soft- 
parton radiation, that need to be summed to all orders to obtain a reliable result. 

In particular, large mass logarithms ln(m|^/m^), which appear in the 6-quark energy 
spectrum, can be resummed using the approach of perturbative fragmentation func- 
tions [8], which expresses the energy spectrum of a heavy quark as the convolution of 
a coefficient function, describing the emission of a massless parton, and a perturbative 
fragmentation function D{mh, fip), associated with the fragmentation of a massless par- 
ton into a massive quark. The method of perturbative fragmentation can be used as long 
as the heavy-quark mass m is much smaller than the hard scale of the process Q, i.e. 
m <^ Q. Given the current limits on the Higgs mass [2], the perturbative fragmentation 
approach can certainly be used in if — > bb, since nib ^ ^i?- 

The dependence of D^rrif,, Hf) on the factorization scale /ip is determined by solving 
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the Dokshitzer-Gribov-Altarelli-Parisi (DGLAP) evolution equations [9,10], once an 
initial condition at a scale /iqf is given. The universality of the initial condition of the 
perturbative fragmentation function, first computed in [8], has been proved in a more 
general way in [11]. 

Moreover, both coefficient function and initial condition of the perturbative fragmen- 
tation function contain terms that become large once the b energy fraction xj, gets close 
to 1. Such terms correspond to soft-gluon emission, and need to be resummed. These 
contributions are process-dependent in the coefficient function and process-independent 
in the initial condition of the perturbative fragmentation function. The resummation of 
soft contributions to the coefficient function oi H ^ bb will be investigated below. 

Finally, in order to describe the 6-quark non-perturbativc fragmentation into b- 
flavoured mesons or baryons B, some phenomenological hadronization models can be 
used. Relying on the universality of the hadronization mechanism, we can tune such 
models to data on B production in e~^e~ annihilation data from LEP or SLD and use 
them to predict the B-hadron spectrum in Higgs decay. Alternatively, we can use exper- 
imental data on the moments of the B spectrum in e^e^ processes, fit the moments of 
the non-perturbative fragmentation function and predict hadron-level moments in Higgs 
decay. 

The plan of this paper is the following. In Section 2, I describe the calculation of 
the next-to-leading order (NLO) MS coefficient function. The approach of perturbative 
fragmentation and the resummation of coUinear logarithms In (m|^/m^) will be discussed 
in Section 3. Section 4 describes the implementation of soft resummation in the coeffi- 
cient function. In Section 5, 1 shall present results on the 6-quark energy spectrum in 
top decay and investigate the effect of soft and coUinear resummation. In Section 6, 
hadron-level results in and N spaces are presented, while Section 7 summarizes the 
main results and gives some concluding remarks. 



2 NLO coefficient function 

I consider Higgs decay into bb pairs at next-to-leading order (NLO) in the strong coupling 
constant as 

H{ph) - b{p,)bipi) {g{pg)) , (1) 

and define the variables: 

The quantities Xf, and Xg are the normalized energy fractions of b and g in the Higgs 
rest frame. As they are expressed in the form of Lorentz-invariant quantities, they can 
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be computed in any frame, provided that the four components of the momenta of b, g 
and H are known. 



In the framework of perturbative fragmentation functions, since mi, <^ mn, one can 
write the differential width for the production of a massive b quark in Higgs decay via 
the convolution: 



]_dri 

Tn dxh 



{xb,mH,mb) 



1 dz 

xt Z 



^ ' - To dz 

+0{{mb/mHr] 



1 dTi 



[z,mH,H,HF) 



MS 
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MS 



Xb 



(3) 



In Eq. 0, dVi/dz is the differential width for the production of a massless parton i in 
Higgs decay with an energy fraction 2, Di{x, fiF,^b) is the perturbative fragmentation 
function for a parton i to fragment into a massive b quark, fi and fip are the renor- 
malization and factorization scales, Tq is the width of the Born process H bb. The 
term O {{mb/rriHy), with p > I, represents contributions which are power-suppressed 
for nib ^H- 

In this section I discuss the computation of the coefficient function (l/ro)(if j/d^;, for 
the production of a massless b in the MS factorization scheme. I shall neglect secondary 
bb production from gluon splitting and limit myself to considering the perturbative 
fragmentation of a massless b into a massive b. In the summation on the right-hand side 
of Eq. (jS)) I shall have only the i = b contribution. 

I regularize ultraviolet, soft and collinear singularities in dimensional regularization, 
and define the parameter e, which is related to the number of dimensions d via c? = 4 — 2e. 
In the computation of dTb/dz, care is to be taken about the treatment of the Yukawa 
coupling of the Hbb vertex. In fact, if the bare Yukawa coupling yb were used, the result 
would be the following: 



dz 




[PU^) - SCpSil - z)] 



+ Giz) 



(4) 



In Eq. (jlj, fir is the regularization scale, remnant of the regularization procedure, G{z) 
is a function, independent of e and fir, whose expression will be detailed later, 7^; = 
0.577 ... is the Euler constant, Cp = 4/3, Pqq{z) is the Altarelli-Parisi splitting function: 



l + z' 
1 - z 



(5) 
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Also, in Eq. (0} I have factorized the four-dimension Born width: ^ 

To = (6) 

Equation (0} shows that a pole 1/e is still present in the differential width. The contri- 
bution proportional to Pqq{z) is associated with colliner radiation from the massless b 
quark and needs to be subtracted to give the coefficient function. Analogous contribu- 
tions have been found, e.g. in the computation of the differential rate of other processes 
such as e~^e~ annihilation [8] or top decay [12]. The additional term, where the pole 1/e 
multiplies the quantity ~ 3as CfS{1 — z), instead has ultraviolet origin, and is charac- 
teristic of Higgs decay and of the scalar nature of the coupling of the Higgs to quarks. 
In fact, unlike the vectorial current, which is conserved, the scalar current is anomalous. 
Because of this extra term, if one naively calculated the total width integrating Eq. (@)), 
one would still find the 1/e pole, which is clearly unphysical. ^ 

In order to get a physical result, one must renormalize the Yukawa coupling. In the 
MS renormalization scheme, the renormalized coupling is related to yb via (see, 

for instance, the discussion in Appendix D of Ref. [13]): 



Vbif^) = yb<l 



(1 TTh^ 
-7 + ^'- + '°«S^^ 



+ . (7) 



One can therefore reabsorb the term ~ 3 asCpSll—z) of Eq. in the MS -renormalized 
coupling yb{iTLH) and evaluate the Born width in © in terms of ybifrLn)- 

In the SM the Yukawa coupling is proportional to the quark mass and, as discussed 
in [13], Eq. ((7j) is consistent with expressing yb{mH) in terms of the MS h mass, mbimu): 

_ grhb{mH) 
yb{mH) 



y/2mw 

where mw is the W mass and g is the coupling constant of SU{2). 

Furthermore, in Refs. [14, 15], where the authors performed the computation with a 
massive h quark, it was shown that a large logarithm ~ a5ln(m^/mf) appears in the 
NLO total rate if one uses the 6-quark pole mass in the coupling. Such a mass logarithm 
can be reabsorbed in the MS mass rhb{mH), which is related to the pole mass rrib by: 



mbiniH = nib 




(9) 



^In order to get the correct finite term in Eq. (@J, one should compute the LO width to 0{e) in 
dimensional regularization. The d-dimension width Td is related to Fq by: Td = ro(47r/m|f)'^[l + e(2 — 
7)]. Equations and following account for the correct finite term. 

'^Pqq{z) being an overall plus distribution, the integral of the term proportional to Pqq(z) over z is 
clearly zero. 
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Subtracting the term ~ Pgq{z){—1/ e + 7 — log47r) from Eq. (jH) and giving an explicit 
expression to the function G{z), one will get the MS H ^ bb coefficient function: 



1 dTh 
f n dz 



MS 



6(1 -z) + 



271 



l + z' 
1 - z 



+ 



-7r^ + -j5(l-z) + l-^-- 
l + z)[ln(l-2)+21n^]+6- 



In^ 

2(1-2;)+ 

In z 



l-z) 



In 2; /ln(l — z) 

2 + 2 ^ 

l-z \ l-z 



(10) 



In Eq. fIlOp I have accounted for the renormalization of the Yukawa coupling and denoted 
by Fq the LO width in terms of ybifnu)- Also, in (fTUI) the factorization scale ^f will have 
to be taken of the order of the Higgs mass, in such a way that the logarithm ln(m^//i|n) 
does not become too large. 

In the following, I shall often make use of the MS coefficient function in Mellin 
moment space Ftv, which is defined by: 



N 



dz z = — —I 
Fn dz 



z . 



In moment space Eq. (fTUI) reads: 

as{^J)CF 



N 



= 1 + 

2 2 



27r 

3 1 

2 ^ iV 



+ 



In 



m 



H 



N + 1 m {N + iy 

[7 + ^o(iV + 2)] 



4^i(iV) + l[7 + V^o(iV + l)] + ^^ 



+ -S,{N + l) + Sl{N -1) + S2{N -1] 



:iii 



(12) 



In Eq. fjl2|) . I have introduced the polygamma functions, ipk{x), which are related to the 
Euler gamma function F(x) through:^ 



rf'^+MogF(x) 
dx'''^^ 



(13) 



■^Reference [12] presents a typing mistake, since ipk{x) is there defined as the fc-th derivative of r(x). 
The numerical results of [12] are nonetheless correct. 
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Equation (fT^ contains also the following combinations: 

Si{N) = ^o(iV + l)-V^o(l), (14) 
S,(N) = -V;i(iV+l)+V^i(l). (15) 

In moment space the convolution Q can then be rewritten as 

TN{mH, rrih) = T^imH, /i, HF)Db,Nil^F, rrih), (16) 

where Tf^^mHjmb) and Db^j^{iip,mb) are the moments of the massive differential rate 
and of the perturbative fragmentation function, respectively. 



Perturbative fragmentation 
and collinear resummation 



The perturbative fragmentation function Db{x, ^Pi^b) introduced in Eq. expresses 
the transition of a massless h into a massive h. Its value at any scale can be obtained 
by solving the DGLAP evolution equations [9, 10], once an initial condition is given. 
As shown in [11], as long as contributions proportional to powers of {uib/mHY can be 
neglected, the initial condition of the perturbative fragmentation function at a scale 
liQF is process- independent. The NLO initial condition in the MS factorization scheme 
reads [8]: 



D'^\xb,fJ,oF,mb) =6{1-Xb) + 



27r 



1-Xb 



fln/if -21n(l 



Xb 



1 



(17) 



The authors of Ref. [16] have recently calculated D™\x, ^oF,^b) to next-to-next-to- 
leading order (NNLO), i.e. up to C(a|). For the purpose of this paper, where the 
coefficient function has been calculated to NLO, the perturbative fragmentation will be 
used to NLO as well. 

The solution of the DGLAP equations in the non-singlet sector, for the evolution 
from the scale fiQp to fip, is given by: 



Db,N{f^F,mb) 



D'^'^{HoF,fnb)exp 



P. 



(0) 



N j^"s(/i 



of) 



2nbn 



asifJ-F 



+ 



«s(/^of) - «5(/^ 



P. 



(1) 



27r6i 



N 



P 



(0) 



N 



In Eq. (fTH|) . D^^^]^{ij,Qp,mb) is the A^-space counterpart of Eq. (fTTj) : P^^ and P^y^ are 
the Mellin transforms of the LO and NLO Altarelli-Parisi splitting functions, and their 
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expression can be found in [8]; 60 and 61 are the first two coefficients of the QCD (3- 
function 

_ 33 - 2nf _ 153 - 19nj 
which enter in the following expression for the strong coupling constant at a scale Q^: 



(n^\ 1 fi fei In [In(gVA^)] 1 

In ()19|) . is the number of active flavours. Equation (|TH)) resums to all orders terms 
containing In (/i|,//iQp). In particular, leading (LL) {agln^lfij^/ flip)) and next-to-leading 
(NLL) (ag ln"'~^(/i|,//iQp)) logarithms are resummed. For an evolution from /iof — 
rrib to — mass logarithms ln(m^/mf) are hence resummed to NLL accuracy 
(collinear resummation). Moreover, setting /ioF — f^ib in Eq- (fT7j) prevents the logarithm 
ln(/iQ^/m^) from getting too large. 

If the calculation were performed with a massive h quark, along the lines of Refs. [14, 
15], contributions ~ a5Pgq(xfe) ln(m|^/m^), equivalent to the collinear pole in massless 
approximation, would be found in the differential spectrum. Hence, using the DGLAP 
evolution equations allows the large mass logarithms appearing in the massive compu- 
tation to be resummed. 

Before closing this section, I wish to point out that, for the purpose of factoriza- 
tion and collinear resummation, using in D™'(xf,, /iof, '^b) and in the DGLAP evolution 
equations the pole or the MS h mass is not as essential as it is in the Yukawa coupling 
((21). In fact, both mass definitions lead to the same results within the given LL or NLL 
logarithmic accuracy. In the following, I shall assume that in Eq. (fTTj) nib is the pole 
mass and will let /ioF run in the range mb/2 < /ioF < 2mb in order to investigate the 
scale dependence of the prediction. ^ The bottom pole mass has also been used in the 
phenomenological analyses of Refs. [11,12,17], within the framework of perturbative 
fragmentation functions. 



4 Soft resummation 



The MS coefficient function (fTIHl and the initial condition of the perturbative fragmen- 
tation function (fTTjl present terms ~ 1/(1 — Xb)+ and ~ [ln(l — Xb)/{1 — Xb)]+ that 
become large once Xb —>■ 1, which corresponds to soft-gluon emission. In moment space, 
they correspond to single ~ In and double logarithms ~ In^ A^ for large values of the 

^If one wanted to use the MS b mass 'fhb{fJ,m), fJ^m should be taken of the order of rrib rather than 
mn- In fact, according to the factorization formula there is no dependence on ttih and on the 
hard-process variables in the perturbative fragmentation function Dt{xb, HPjiTib)- 
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Mellin variable N. Such contributions are process-independent in the initial condition 
of the perturbative fragmentation function, and have been resummed in [11] in the NLL 
approximation . 

In this section I would like to present the results for soft resummation in the MS 
coefficient function. First, it is instructive to write Eq. (|12j) for large A^: 



1 + 



In^ + - + 27 



21n^ 



In AT 



N 



where K{mH, fJ'p) contains terms that are constant with respect to A^: 



K{mH,f^F) 



In ^ 



3 
2 



-lE + 1e- 



(21) 



(22) 



Furthermore, to get Eq. ()2H). I have used the large- A^ expansions of the polygamma 
functions: 



^i(iV) 
Si{N) 

S2{N) 



InA^ + O 



In AT + 7s + 
/ 1 

— + C — 

6 VA^ 



(23) 
(24) 
(25) 
(26) 



LL and NLL soft contributions to the MS coefficient function can be resummed following 
standard methods [18,19]. 

In particular, all the steps that led to NLL resummation in the coefficient function 
for e+e~ qq processes in Ref. [11] can be repeated, thus obtaining the resummed 
coefficient function for H bb. In fact, in both processes one has in the final state 
two massless partons, as the b and b are in the coefficient function, which are able to 
radiate soft- as well as coUinear-enhanced radiation. The coefficients of In A^^ and In A^ 
in Eq. (^1]) are indeed the same as in the large-A^ expansion of the e"*"e~ coefficient 
function [11]. 

The resummed coefficient function can be obtained from the e~^e~ one, replacing the 
centre-of-mass energy squared with m'jj. One obtains: 



In A 



N 



dz- 



1 



1 - z 



■A 



2 



as 



(27) 
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In Eq. fjTfjl . the two integration variables are z = \ — Xg and I? = {ph + ^9)^(1 ~ z), 
as in [18]. In soft approximation, z ~ Xb'-, for small-angle radiation, k'^ ~ k'^, the gluon 
transverse momentum with respect to the 6-quark direction. 

The function B{as) is associated with the radiation emitted by the unobserved 
massless parton, namely the h if one observes the h. The argument of B[as) is the 
invariant mass of the unobserved jet, i.e. [pi + PgY — ^'h{^ ~ ^) soft approximation. 

Moreover, Eq. (jTfjl is formally equal to the resummed MS coefficient function in 
Drell-Yan and deep inelastic scattering (DIS) with light quarks [20]. For processes with 
massive quarks, such as top quark decay [17] or heavy quark production in DIS [21], 
recently provided with NLL soft resummation, the function B{as)/2 should be replaced 
by a different one, called S{as) in Refs. [17] and [21], which is characteristic of processes 
with heavy quarks and expresses soft radiation, which is not coUinear-enhanced. The 
function A{as) in Eq. (|77|) can be expanded as a series in as as: 



f28) 



The first two coefficients are mandatory to resum the coefficient function to NLL accu- 
racy [18]: 

A^^^ = Cf, (29) 

.2\ cr 



= -Cp 

2 



Q7 TT' 
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where Ca = 3. Likewise, the function B{as) can be expanded: 

and, to NLL level, only the first term of the expansion is kept: 

The integral over z can be performed making use of the replacement [18]: 

1^-0(1 



-lE 



N 



(30) 



(31) 



(32) 



(33) 



where is the Heaviside step function. The function Ajy can be expressed in the usual 
form: 

Ajv(m^^,^,/i^) =exp[lnAr^«(A)+^(2)^^^^^^^)] ^ (34) 

with 

A = 6ott5(/i)lniV. (35) 
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In Eq. flH4|l the term In Ng^^\\) accounts for the resummation of the leading logarithms 
a^\n^~^^ N in the Sudakov exponent, and the function g'-'^\X, hf) resums NLL terms 
In" N. Functions g^^^ and g^'^^ can be obtained by simply setting Q"^ = rn?^ in Eqs. (34) 
and (35) of Ref. [11]; the result is not here reported for the sake of brevity. 

Following [11,12,21], I include in the Sudakov-resummed coefficient function the 
constant terms K{mH, fJ^p) defined in Eq. and obtain: 



X exp 



27r 

lniV^«(A) + (7(2)(A,^,;,^) 



(36) 



Finally, the resummed result is matched to the fixed-order one, which will yield a pre- 
diction valid at Xfe < 1 as well. The NLO result is added to Eq. (jHUj) and, in order to 
avoid double counting, the 0{as) term of the resummed result is subtracted: 



+ 



f^(mH,/i,yUF) - 
f Ar(m/^,/i,/ii7) 



as 



as 



(37) 



where [f ^flas [^N]as respectively the expansion of Eq. ^T7\ up to 0{as) and 
the full fixed-order coefficient function at 0{as) (fT^. 



5 Bottom-quark energy distribution 

In this section I present the energy spectrum of bottom quarks in Higgs decay, accord- 
ing to the calculation above described. I shall investigate the phenomenological effect 
of collinear and soft resummation, and the dependence of the prediction on factoriza- 
tion and renormalization scales. The results, which have been presented in analytic 
form in Mellin space, will be inverted numerically to Xb space choosing the integration 
contour according to the minimal prescription [22]. I shall plot the normalized rate 
{l/r){dr /dxb), where F is the NLO H ^ bb width, calculated in Refs. [14,15]. I shall 
set niH = 120 GeV, = 5 GeV, ri/ = 5, A = 200 MeV. 

In Fig. ^ I have plotted the energy spectrum of the 6-quark, according to the NLO 
massive calculation (dashed line), which can be obtained from Eq. (jS)) without evolving 
the perturbative fragmentation function, including collinear resummation (dotted line) 
and with both collinear and soft resummations (solid). I have set = = ^h, 

We can notice a remarkable effect of both collinear and soft resummations. The 
fixed-order calculation lies below the resummed ones and grows as ^ 1 because of 



10 



a behaviour ~ 1/(1 — Xb)- The colhnear-resummed spectrum exhibits instead a sharp 
peak at large Xb, after both soft and coUinear resummations, the distribution is further 
smoothed and presents the Sudakov peak at Xb — 0.97. 

Although one implements both coUinear and soft resummations, the prediction is 
still not reliable at very small and large Xb- In fact, even the resummed distributions 
become negative once Xb approaches or 1. At small Xb, the coefficient function contains 
a term ~ logx;, which has not been resummed yet; at very large Xb, one starts to get 
sensitive to missing non-perturbative power corrections. The range of reliability of a 
purely perturbative computation is typically x;, < 1 — A/rrib [11]. 

In Figs. I2E I show the dependence of the prediction on the factorization and renor- 
malization scales that appear in the calculation, i.e. fi and /zp in Eq. (|TO|l . /xq and /zof 
in (|T7jl . I let fi and fip assume the values mH/2, mu and 2m h] as for the scales /iq 
and yUoF, as discussed in Section 3, they will be taken equal to nib and 2mb. The 

dependence on the scales is logarithmic, hence we expect it to be more visible once the 
scales vary around small values, i.e. around irib rather than mH- 

CoUinear resummation is turned on in all figures, but plots are with or without soft- 
gluon resummation. The general feature of these results is that after both coUinear- 
and soft-enhanced terms are resummed, the final prediction exhibits very little scale 
dependence. In particular, we note that the role of soft resummation is crucial to 
weaken the dependence on scales iip (Fig- 12)), /^o (Fig- El) and especially /xof (Fig. E))- 
The dependence on the renormalization scale is instead very weak even resumming 
only coUinear logarithms, as shown in Fig. |3J A milder sensitivity to renormalization 
and factorization scales corresponds to a reduction of the theoretical uncertainty on the 
result and is a meaningful effect of collinear and soft resummations. 

Moreover, I would like to compare the 6-quark energy distribution in three processes 
that have been investigated within the framework of perturbative fragmentation func- 
tions and provided with NLL collinear and soft resummations, namely e^e~ — >■ hh [11], 
top decay t bW [12, 17] and, in this paper, Higgs decay. Because of the universality 
of the fragmentation process, possible differences in such spectra will be related to the 
different coefficient functions and mass scales involved. Furthermore, I shall consider 
e~^e~ at a/s = 91.2 GeV, the centre-of-mass energy of LEP I and SLD, and, in order 
to test the effect of the coupling, vectorial in e~^e~ hh and scalar in H ^ hh, also 
at ^/s = 120 GeV, the default Higgs-mass value throughout this paper. In top decay 
I shall set rrit = 175 GeV and mw = 80 GeV. From Fig. IHlone learns that the shapes 
of the three distributions exhibit some differences. The h spectrum in H decay is the 
highest at small Xb and the lowest at large Xb, in top decay it is shifted toward large Xb 
and peaked very close to 1. The e~^e~ — >■ hh prediction lies within the other two at small 
and very large Xb- The effect of different values of y/s in the e+e~ process is visible 
mainly around the Sudakov peak; setting y/s = uih makes the spectrum more similar 
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X 
t3 




Figure 1: 6-quark energy distribution in Higgs decay according to the NLO massive 
calculation (dashed line), including NLL coUinear resummation (dots) and both NLL 
coUinear and soft resummations (solid). I have set: fj. = fip — ttih, A'-o = A*of = iTib- In 
the inset figure, the same curves are shown at large Xf, and on a logarithmic scale. 



XI 

73 




Figure 2: Dependence of the b spectrum on the factorization scale iip, with (solid fines) 
and without (dots) NLL soft resummation. The other scales are fixed at /i — mn and 
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Figure 3: As in Fig. |21 but for different values of the factorization scale /^of- The other 
scales are fixed ai = = uih and fiQ = nib. 




Figure 4: As in Fig. |21 but varying the renormalization scale /.i. The other scales are 
fixed a.t fip = rriH and /iof = /^o = f^h- 
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Figure 5: As in Fig. but for different values of the renormalization scale /xq. The other 
scales are fixed dX n = hf = and fioF = ^b- 

to the Higgs decay one. 

Before closing this section, it is interesting to investigate the dependence of the 
resummed prediction on the Higgs mass. In Fig. [71 we plot the NLL collinear- and soft- 
resummed spectrum for rriH = HO, 120 and 130 GeV. We note that the spectra change 
very little with the chosen value of mn, and some small effect is only visible around the 
Sudakov peak. 



6 Hadron-level results 



I shall now present results on the energy spectrum of 6-flavoured hadrons B in Higgs 
decay. Similarly to the parton-level analysis, for a hadron of momentum ps one defines 
the B normalized energy fraction in the Higgs rest frame: 

Xb = ?^E^_pL_ (38) 



m 



H 



Up to power corrections, one writes the hadron-level spectrum as the convolution of the 
parton-level one, calculated above, with the non-perturbative fragmentation function 
D"'P{x), associated with the hadronization of the b quark into the B hadron: 



r dx 



XB,mH,mb) = - [ ——^{z,mH,mb)D''f(—\ (39) 
B r JxB z dz \ z J 
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Figure 6: 6-quark energy distribution in Higgs decay (solid), top decay (dots) and 
e+e~ bb processes ai ^/s — 91.2 GeV (dashes) and 120 GeV (dot-dashes). All 
predictions are given by a NLO calculation provided with NLL coUinear and soft resum- 
mations. In Higgs decay I have set n = hf = rriH, in top decay ^ = — rrit — 175 
GeV, in e'^e~ annihilation ii = fip = AH plots are for fiQ = fiQp = ^h- 



X 
t3 




Figure 7: 6-quark energy distribution in Higgs decay, including NLL coUinear and soft 
resummations, for a Higgs mass of 110 GeV (dashes), 120 GeV (sohd) and 130 GeV 
(dots). 
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The non-perturbative fragmentation function is to be extracted from experimental data. 
As the Higgs particle has not been discovered yet, one obviously does not have any 
data on B production in Higgs decay. However, relying on the universality of the 
hadronization transition, one can fit some hadronization models to data on B production 
in e+e~ experiments and use them to predict the B spectrum in H ^ bb processes. In 
order for the fitting procedure to be consistent, the e^e~ bb perturbative process is 
to be described as done for Higgs decay. I shall have to use NLO coefficient functions, 
NLL collinear and soft resummation, and consistent values for the scales appearing in 
the calculation, e.g. /z = ^/s if I set /i = mn in H ^ bb. In fact, it was recently found 
[23, 24] that such a consistency is crucial if one wishes to make use of non-perturbative 
information taken from e~^e~ data to describe the data on 5-hadron production at 
hadron colliders, such as the Tevatron. 

In this paper, as for the non-perturbative fragmentation function, I shall consider 
the following three hadronization models: a power law with two parameters 

the model of Kartvelishvili et al. [25] 

= (l + 5)(2 + 5)(l-x)x^ (41) 

and the Peterson model [26] 

^"^(^; = ^ -. (42) 

x[l — l/x — e/(l — x)Y 

In Eq. fl40j) . B{x, y) is the Euler beta function; in ()42|) A is a normalization constant. The 
parameters a, /9, 5 and e are to be extracted from experimental data. In Ref. [17], such 
models have been fitted to ALEPH [27] data on B mesons and it was found that models 
()40|) and ()41|) yield very good fits, while the Peterson model is marginally consistent. In 
Ref. [28] , also the SLD data [29] were considered and it was found that Eqs. (jlUj) and (PT|) 
lead to good fits, and Eq. is instead unable to reproduce the SLD data. Moreover, 
using a soft-resummed perturbative calculation was essential to describe the SLD data. 
However, discrepancies were found in Ref. [28] between the best-fit parameters a, (3 and 
5, according to whether one fits the models to ALEPH or SLD. Unlike ALEPH, the SLD 
data contain some 6-flavoured baryons, mainly the A;,, but it is a pretty small fraction of 
the whole sample; hence, it may not be correct to conclude that the differences reported 
in [28] are due to the baryons reconstructed in SLD. Moreover, detailed analyses and 
comparisons on ^-fragmentation in the four LEP experiments and SLD are currently 
missing. 

Here I shall try instead a combined fit of both ALEPH and SLD data samples and 
investigate whether one is able to find a suitable parametrization of the hadronization 
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models (jin|) - (jl^ that would be consistent with both experiments. As in Refs. [12,17,28], 
I shall consider data for 0.18 < < 0.94, so as to avoid data points close to or 1, 
where the presented calculation is not reliable. Furthermore, when doing the fits, I shall 
neglect the correlations among data points and sum statistical and systematic errors in 
quadrature. 



a 


0.90 ±0.15 




16.23 ± 1.37 


x'(a,/3)/dof 


33.42/31 




17.07 ±0.39 


x'(^)/dof 


33.80/32 


e 


(1.71 ±0.09) X 10-3 


x'(e)/dof 


166.36/32 



Table 1: Results of combined fits to e^e' bb data from ALEPH and SLD col- 
laborations, using NLO coefficient functions, NLL DGLAP evolution and NLL soft- 
gluon resummation. I have set A = 200 MeV, fiQp = fio = rrib = 5 GeV and 
fip = fJ' = \/s = 91.2 GeV. a and /3 are the parameters in the power law (HOJ, 6 
refers to (gH), e to (gH). 

In Table Q the best-fit parameters are quoted, along with the P^r degree of free- 
dom. Both power law and Kartvelishvili models fit the data quite well, while the Pe- 
terson non-perturbative model is unable to describe them. It is also interesting to 
investigate whether the implementation of soft resummation has an impact on the fit. 
Table 121 shows the best-fit parameters obtained without resumming soft-gluon contribu- 
tions to the perturbative calculation of e'^e" bb. It should be noted that the fit gets 
worse and no hadronization model is capable of yielding a reasonably small x^/dof. 

In Fig.lHlI present a prediction on the xb distribution in Higgs decay, using the NLL 
coUinear- and soft-resummed perturbative calculation, and the best-fit parameters of 



a 


1.01 ±0.12 


/3 


14.73 ± 0.94 




64.19/31 


6 


14.67 ±0.30 


x'(5)/dof 


64.20/32 


e 


(2.76 ±0.14) X 10-3 


x\e)/doi 


287.55/32 



Table 2: As in Table CI but without soft resummation in the parton-level calculation. 
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0.2 0.4 0.6 0.8 1 

Figure 8: i?-hadron spectrum in Higgs decay, modelling the hadronization according 
to the power law ()4()j) (solid) and the Kartvelishvili model (dashes). Plotted are 
the edges of bands at one-standard-deviation confidence level for the non-perturbative 
parameters a, (3 and 5, as reported in Table ^ In the perturbative calculation, the 
scales have been set to: ^ = fip = t^H = 120 GeV, /io = /ioF = = 5 GeV. 

models pn|) and PT|) quoted in Table [H I discard the Peterson model as it does not 
acceptably describe the considered data sample. In order to account for the errors on 
the best-fit parameters given in Tabled for each model I plot a band corresponding to 
a prediction at one-standard-deviation confidence level for the fitted non-perturbative 
parameters a, j3 and S. 

Fig.lSlshows that the two predictions yielded by models (jlUj) and (PT|) are statistically 
consistent. In fact, from Tabled one learns that, within the error range, the parameter 
a of Eq. (^njl is consistent with 1 and /? of Eq. pUjl is consistent with 6 of Eq. pT|l . It 
is therefore reasonable that the predictions of the power law with two parameters and 
of the Kartvelishvili model agree. 

As I did for the analysis at parton level, I present in Fig. |H1 the S-hadron spectra in 
three different processes, i.e. e~^e~ annihilation, Higgs and top decay, using everywhere 
a NLO and NLL resummed perturbative calculation and the hadronization model ()4U|). 
with the best-fit parameters reported in Tabled As in the parton-level analysis, I make 
consistent choices for the scales involved and consider two centre-of-mass energies for 
the electron-positron process, i.e. ^/s = 91.2 and 120 GeV. For each process I plot a 
band corresponding to a prediction at one-standard-deviation confidence level for the 
parameters a and /? of Eq. . 
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Figure 9: 5-hadron spectra in Higgs decay (solid), in top decay (dashes) and in e~^e~ 
annihilation at ^/s = 91.2 GeV (dotted) and 120 GeV (dot-dashed), according to the 
hadronization model (j4Uj) . at one-standard-deviation confidence level for the parameters 
a and (3. Factorization and renormalization scales are chosen as in Fig. IHl 

Figure IHl shows that the three spectra are statistically different and that the result 
is pretty similar to the one already found at parton-level in Fig. which is reasonable 
as we are convoluting the Xb distribution with the same non-perturbative fragmentation 
function. Setting -^i = rriH makes the e"'"e~ bb spectrum closer to the if-decay 
one, especially at middle and large values of xb, though meaningful differences are still 
present at small xb and around the peak. 

Before closing this section, I would like to present results in moment space, using 
the experimental moments on B production in e^e~ annihilation from the DELPHI 
Collaboration [30], along the lines of [17]. The advantage of working in moment space, 
as suggested in [23], is that one does not need to rely on any specific functional form of 
the non-perturbative fragmentation function. As done in [17], I have quoted in Table El 
the first four moments on 5-production at DELPHI, the computed moments of e^e~ and 
Higgs decay perturbative calculations, the extracted moments of the non-perturbative 
fragmentation function and the predictions for the moments of the B spectra in H bb. 

In moment space convolutions are turned into ordinary products, and we therefore 
have: a§ = a^D^, = V'^^D^ = Vfja^/a^^, where and Ftv are the moments 
of e"'"e~-annihilation cross section and Higgs-decay width at hadron or parton level. 
Comparing Table El with the results of [17], it should be noted that all considered 
moments in H decay are lower than those in e+e~ processes and in top decay. This 
result is consistent with the spectra in x^-space presented in Fig. El 
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(x) 










0.7L33±0.(J(J52 


0.5 101±0. 0061 


0.1236±0.0065 


0.3 10(J±0. 0064 


e+e- NLL 


0.7801 


0.6436 


0.5479 


0.4755 




0.9169 


0.8392 


0.7731 


0.7163 


H-decay NLL L^ 


0.7578 


0.6162 


0.5193 


0.4473 


ff-decay 


0.6948 


0.5171 


0.4015 


0.3204 



Table 3: Experimental data for the moments from DELPHI [30], the resummed e"'"e~ 
perturbative calculations for ct^vt [H]) extracted non-perturbative contribution D^^\ Using 
the resummed perturbative result F^, a prediction for the moments in H ^ bb processes 
is given. The experimental error should be propagated to the final prediction. 

7 Conclusions 

I have considered bottom-quark fragmentation in Standard Model Higgs decay H ^ bb 
within the approach of perturbative fragmentation functions. I have computed the 
MS NLO coefficient function and resummed collinear logarithms ln(mf^/m^) by using 
the DGLAP evolution equations in the NLL approximation. In the calculation of the 
coefficient function, the use of the MS -renormalized Yukawa coupling turned out to be 
essential. 

I have resummed NLL soft terms in the coefficient function and matched the re- 
summed result with the exact NLO one. Soft resummation in the coefficient function 
has been combined with the process-independent NLL soft resummation in the initial 
condition of the perturbative fragmentation function. 

I have presented the 6-quark energy spectrum in Higgs decay, which exhibits a re- 
markable effect of the implemented collinear and soft resummations. In particular, soft 
resummation smoothens the distribution at large and the dependence of the prediction 
on factorization and renormalization scales turns out to be very little. The dependence 
of the predicted spectra on the Higgs mass is very small as well. I have also compared 
the 6-energy spectrum in H ^ bb with the ones yielded by top decay and e'^e~ — > bb, 
which have been provided with NLL collinear and soft resummations in the framework 
of perturbative fragmentation. 

I have then considered 6- flavoured S-hadron production in Higgs decay in both xb 
and moment space. I have fitted a few hadronization models to ALEPH and SLD data 
in x^-space and used the best-fit parameters to predict the B spectrum in H decays. For 
this procedure to be consistent, the perturbative processes e"'"e~ — > bb and H ^ bb have 
been described using the same kind of calculation. The fits have shown that the power 
law with two tunable parameters and the Kartvelishvili model give good combined fits 
of ALEPH and SLD data, while the Peterson model is not capable of reproducing the 
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data. Moreover, implementing soft-gluon resummation in the coefficient function has 
been essential to be able to describe the data. The 5-hadron spectra in Higgs decay 
according to the power law and to the Kartvelishvili model are in statistical agreement. 
The comparison of B spectra in H ^ bb, t bW and e~^e~ — > bb exhibits similar 
features to the parton-level one. I have also made predictions for the ffist four moments 
of the B spectrum in H ^ bb processes, using DELPHI data on the B moments in e^e~ 
annihilation. 

In conclusion, the presented calculation allows the performance of precise predictions 
for 6-quark and i?-hadron production in Higgs decay and could be applied for analyses of 
Higgs phenomenology at the Tevatron and ultimately at the LHC. The computed NLO 
matrix elements, along with NLL DGLAP evolution and soft resummation, can be used 
for further investigations of other observables in Higgs decay. In particular, studies of the 
transverse momentum distributions of b quarks and 6-flavoured hadrons are under way. 
Moreover, as the NNLO initial condition of the perturbative fragmentation function 
has recently been computed [16], the coefficient function may be calculated to NNLO as 
well and a resummation of collinear and soft terms to next-to-next-to-leading logarithmic 
accuracy can be made. 
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